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ABSTRACT 

We study secular stability against a quasi-radial oscillation for rigidly rotating stars with soft equations 
of state in general relativity. The polytropic equations of state with polytropic index n between 3 and 
3.05 are adopted for modeling the rotating stars. The stability is determined in terms of the turning-point 
method. It is found that (i) for n ^ 3.04, all the rigidly rotating stars are unstable against the quasi- 
radial oscillation and (ii) for n J; 3.01, the nondimensional angular momentum parameter q = cJ /GM^ 
(where J, M, G, and c denote the angular momentum, the gravitational mass, the gravitational constant, 
and the speed of light, respectively) for all marginally stable rotating stars is larger than unity. A semi- 
analytic calculation is also performed, and good agreement with the numerical results is confirmed. The 
final outcome after axisymmetric gravitational collapse of rigidly rotating and marginally stable massive 
stars with g > 1 is predicted, assuming that the rest-mass distribution as a function of the specific 
angular momentum is preserved and that the pressure never halt the collapse. It is found that even for 
1 < g ^ 2.5, a black hole may be formed as a result of the collapse, but for q 2.5, the significant 
angular momentum will prevent the direct formation of a black hole. 

Subject headings: black hole physics - relativity - hydrodynamics - stars: rotation 



1. INTRODUCTION 

At the final stage of their evolution, massive stars with 
initial mass 8Mq form an iron core. Because of high 
temperature and density, the photodissociation of the iron 
to lighter elements proceeds rapidly (Shapiro & Teukol- 
sky 1983). In the photodissociation, the thermal pressure 
is consumed, decreasing the adiabatic index F below 4/3. 
Consequently, the core is destabilized, resulting in collapse 
to a neutron star or a black hole. 

Recent numerical simulations of the collapse of primor- 
dial molecular clouds have suggested that the first gen- 
eration of stars contains many massive members of mass 
- 100-IOOOMq (Bromm et al. 1999; Abel et al. 2000; 
Nakamura & Umemura 2001). Such massive stars form 
large helium cores that reach carbon ignition. It is known 
that after helium burning, cores of such stars encounter the 
electron-positron pair creation instability, bringing F be- 
low ~ 4/3 and triggering the collapse (Barkat et al. 1967; 
Bond et al. 1984; Kippenhahn & Weighrt 1990; Fryer et al. 
2001). For sufficiently massive stars, of mass M ^ 260Mq, 
such unstable cores form a black hole, while for less mas- 
sive cases, the outcome is a pair-unstable supernova (Fryer 
et al. 2001). 

These examples illustrate that the gravitational col- 
lapses set in when the stability of progenitors changes, 
and hence, a star marginally stable against quasi-radial 
oscillations may be regarded as an initial condition of the 
collapse. Massive stars in nature are rapidly rotating in 
general (Bond et al. 1984). Thus, determining the quasi- 
radial stability for rotating stars is an important subject 
in astrophysics for elucidating a plausible initial condition 
for the stellar collapse. 

In many cases, one simply assumes that progenitors of 
neutron stars or black holes start collapsing when their 
adiabatic index F decreases below a well-known critical 



value, 4/3, against a quasi-radial instability. The number 
4/3 is the correct value only when the progenitor is not 
rotating. If the progenitor is rotating, it is not the valid 
criterion, since angular momentum has the effect of stabi- 
lizing the star against gravitational collapse. In particular, 
the number 4/3 would not be a good indicator of the insta- 
bility for pair-unstable collapse, since the progenitor stars 
are likely to be rapidly rotating (Bond et al. 1984, Fryer 
et al 2001). 

According to a Newtonian theory, slowly rotating stars 
that satisfy the following relation are unstable against a 
quasi-radial oscillation (Tassoul 1978; see also §2): 



3F - 4 - 2(3F - 5)/3 < 0. 



Here, 
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T 



(1) 



(2) 



and T and W are rotational kinetic and gravitational po- 
tential energies, respectively. In the following, we restrict 
our attention to the rigidly rotating case for simplicity. 
Based on our numerical analysis in Newtonian gravity, 
the maximum value of (3 for rigidly rotating stars with 
polytropic equations of state, which is achieved when the 
velocity at the equatorial surface is equal to the Kepler ve- 
locity (i.e., at the mass-shedding limits), is approximately 
written as 



/9> 



n,ax - 0.00902 + 0.124(F-- 



(3) 



Combining equations (1) and (3), we find that for F ^ 
1.328, all the rigidly rotating stars are unstable against 
a quasi-radial oscillation in Newtonian gravity. However, 
even for 1.328 ^ F < 4/3, the stars can be stabilized by 
the effect of rotation. 
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General relativistic gravity destabilizes the rapidly ro- 
tating stars even in the case of F > 1.328. One of the 
purposes of this paper is to numerically determine the cri- 
teria for the onset of the instability. In particular, we 
focus on the secular stability of rigidly rotating stars with 
soft equations of state in equilibrium against a quasi-radial 
oscillation. The stability is determined using the turning- 
point method (Friedman et al. 1988). To model the rotat- 
ing stars, we adopt the polytropic equations of state with 
n ^ 3 (F ^ 4/3). 

From a general relativistic point of view, it is interesting 
to ask how large the nondimensional angular momentum 
parameter q = J/M^ can be for the marginally stable stars 
that may be plausible progenitors of gravitational collapse. 
Axisymmetric hydrodynamic simulations in general rela- 
tivity (Nakamura 1981; Stark & Piran 1985; Piran & Stark 
1986; Nakamura et al. 1987; Shibata 2000) have provided 
many numerical results that show that for g > 1, a black 
hole is not formed after the collapse. That is, the value of 
q may be a good indicator in predicting the final fate of 
the stellar collapse of massive objects. As Baumgarte and 
Shapiro (1999) indicate, the value of q for a marginally 
stable star with F = 4/3 is very close to unity. Cook et 
al. (1994) show that the critical value of q increases with 
decreasing of F for 4/3 < F <, 5/3. This implies that for 
F < 4/3, the critical value is likely to be larger than unity, 
even for the case of 4/3 — F -C 1. The second purpose of 
this paper is to confirm this fact. 

The third purpose is to predict the outcomes of gravi- 
tational collapse of the marginally stable stars. As men- 
tioned above and confirmed below, the value of q for the 
marginally stable stars with F < 4/3 can be larger than 
unity. This suggests that a rotating star with F < 4/3 of 
mass large enough collapse to a black hole in the nonro- 
tating case may not form a black hole because of its sig- 
nificant effect of the angular momentum. In the analysis, 
we assume that (i) the collapse proceeds in an axisymmet- 
ric manner, (ii) the rest-mass distribution as a function 
of the specific angular momentum is preserved, and (iii) 
the pressure never halt the collapse. We indicate that the 
criterion for no black hole formation, q > I, which has 
been suggested so far, is not very good for rigidly rotat- 
ing stars with the soft equations of state. We predict that 
even for 1 < g 2.5, a black hole may be formed as a 
result of gravitational collapse of the rigidly rotating and 
marginally stable stars. However, for q > 2.5, the signifi- 
cant angular momentum will prevent the collapsing stars 
from directly forming a black hole. 

The paper is organized as follows. In §2, we present 
a semi-analytic calculation for determining the stability 
of rotating stars against a quasi-radial oscillation. In §3, 
the secular stability of rotating stars in general relativ- 
ity is numerically determined. In §4, we predict the final 
outcomes of the stellar collapse for very massive stars, as- 
suming that the initial condition is a marginally stable 
star, as determined in §3. Section 5 is devoted to a sum- 
mary. Throughout this paper, the pressure is given by the 
polytropic equation of state as 

P = Kp^, T = l + -, (4) 
n 

where P is the pressure, p the baryon rest-mass density, K 
the polytropic constant, and n the polytropic index, which 



is chosen in the range between 3 and 3.05. We adopt the 
geometrical units G = c = 1, where G and c denote the 
gravitational constant and the speed of light, respectively. 

2. SEMI-ANALYTIC EXPLORATION FOR STABILITY 

To semi-analytically determine the marginally stable ro- 
tating stars with the polytropic equations of state with 
n J> 3, we follow the method described by Zel'dovich & 
Novikov (1971) and Shapiro & Teukolsky (1983) and write 
the total energy as the sum of the internal energy, the 
Newtonian potential energy, the rotational kinetic energy 
(in Newtonian order), and a post-Newtonian correction. 
These terms can be written 

E = kxKMpl/'' - k2M^/^pl'^ 

+k^J^M-^/^pll^ - kiM^^pl/^ (5) 

where k\, k2, k^, and are structure constants that 
are constructed from Lanc-Emdcn functions (Shapiro & 
Teukolsky 1983) and listed in Table 1. M is the total 
mass and Pc the central density. 

To derive equation (5), we assume that the stars are 
spherical and rigidly rotating. In the treatment by 
Zel'dovich & Novikov (1971) and Shapiro & Teukolsky 
(1983), the authors take into account a nonsphericity, as- 
suming that a spherical surface of constant density in the 
nonrotating case transforms into a spheroidal surface en- 
closing the same volume. This approximate method is 
likely to be good for stars with very stiff equations of state 
with n ~ 0. However, this is not the case for the stars with 
n ^ 3, since the density profile is not as uniform as that 
for n = 0. Actually, in their approximate treatment, the 
value of (3 for given values of pc and M is significantly over- 
estimated for n ^ 3. In the case of the soft equations of 
state, the shape of the central region, around which most of 
the mass is concentrated, is nearly spherical even at mass- 
shedding limits. Therefore, we neglect the nonsphericity 
here. In this treatment, the value of J is considered an 
input quantity. Note also that the second post-Newtonian 
term is omitted in the present treatment, in contrast to 
those by Zel'dovich & Novikov (1971) and Baumgarte & 
Shapiro (1999). This is because the compactness of sta- 
ble stars with a soft equation of state is very small, and 
hence, the second post-Newtonian correction produces a 
very small effect on the equilibrium and stability, as illus- 
trated by Baumgarte & Shapiro (1999) for n = 3. 

Taking the first derivative of equation (5) with respect 
to the central density yields a condition for equilibrium, 

= ^ = -k^KMx^^-'^y^ - k^M^'"" 
ox n 

-K2fc5 J^M-s/^x - 2k4M'^/^x, (6) 

where x = fhJ^ . For a stable equilibrium, the second 
derivative of equation (5) has to be positive. Therefore, a 
root of the second derivative marks the onset of a quasi- 
radial instability: 

= ^ = ^ / 3 _ A ^^^^^(3-2„)/n 

ox^ n\n ) 

^Ik^j'^M-^l'^ - Ik^M^I^. (7) 
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By combining equations (6) and (7), we derive a relation 
between M and x for the marginally stable stars as 



M 



3(2n - 3)ki ,^^(3_„)/„ 



-1 3/2 



(8) 



For n = 3, MK~^/'^ is a constant as 4.5548 irrespective of 
X (Baumgartc & Shapiro 1999). For n > 3, M decreases 
with increasing pc for a fixed value of K, and in the limit 
Pc ^ 0, M ^ oo. 

Prom equations (7) and (8), J/M"^ for the marginally 
stable stars (defined as ^mar) can be written as a function 
of the density: 



2 _ ^4 



n^(n-3)fci 
6(2n - 2,fkik5 



(9) 



For r?, = 3, (jmar is a constant as 0.873 and for n > 3, 
(Zmar dccrcascs with increasing pc- In the limit pc — > 0, 
9mar oo. As we show in §3, the maximum value 
of pc along the sequence of the marginally stable stars 
is reached at mass-shedding limits. At that point, g,„ar 
reaches its minimum (<7min)- Thus, the value of Qmar for 
the marginally stable stars is always larger than q^m. for 
n > 3. 

Using equation (6), equation (7) can be rewritten as 
= 3r - 4 - 2(3r - 5) f %^ - ^\ M^'^pI'\ (10) 

V K2 k2 J 

In the present approach, the second term is related to 

/3(= T/W) by 

k2 

Using this relation, equation (10) is rewritten as 

= 3r - 4 - 2(3r - 5) (^/? - ^M^/^y ) . (12) 

The term M'^l'^p'J^ is related to the compactness parame- 
ter M/Rg, where Rg is the radius of the spherical star, as 
follows: 

^=aMy'pl'\ (13) 
tig 

Here, the values of a arc determined by the Lanc-Emden 
function and listed in Table 1. Using this relation, equa- 
tion (12) is written as 



= 3r - 4 - 2(3r - 5)/3 - a' 



M 

Rg ' 



(14) 



where a' = 2(5 — 3r)fc4/(aA:2) which is in the range be- 
tween 6.747 and 7.137 for F = 4/3-1.328. Equation (14) 
is a familiar formula often presented in the standard text- 
books such as Tassoul (1978). For the case in which f3 — 0, 
the equation agrees with the formula derived by Chan- 
drasekhar (1964) for n = 3. 

3. NUMERICAL ANALYSIS FOR STABILITY 

In §2, we determined sequences of the marginally stable 
stars for n ^ 3 (F ^ 4/3) in terms of a semi-analytic cal- 
culation. As shown in this section, they indeed give good 



approximate sequences for the marginally stable rotating 
stars. However, they are not exact sequences after all. 
To determine them precisely, numerical computations are 
necessary. 

Another drawback in the semi-analytic calculation is 
that the maximum value of the angular velocity for a given 
value of J is not determined precisely. For rigidly rotating 
stars, the velocity should be smaller than the Kepler ve- 
locity at the equatorial surface. This restricts the allowed 
region for the rigidly rotating stars: A sequence of stars at 
mass-shedding limits divides the regions where the rigidly 
rotating stars exist and where they do not. To determine 
the sequence at the mass-shedding limits, the configura- 
tion of the stars has to be accurately computed. Thus, 
numerical computation is inevitable. 

3.1. Basic equations 

To study the secular stability of rotating stars against 
a quasi-radial oscillation, the equilibrium solutions in gen- 
eral relativity are computed. We write the energy mo- 
mentum tensor of the Einstein equation for an ideal fluid 
as 



T'"' = phuV + Pg'"' , 



(15) 



where u^^ is the four velocity, h = l+ e + P/pis the en- 
thalpy, £ is the specific internal energy, and g^" is the 
spacctime metric. As in §2, the polytropic equations of 
state are adopted. Using the first law of the thermody- 
namics, the specific internal energy e in the polytropic 
equations of state is written as 



nP 



£ = 



(16) 



As in §2, we pay attention only to rigidly rotating stars 
setting the angular velocity = u'^/u* as a constant. 

With the polytropic equation of state, physical units en- 
ter the problem only through the polytropic constant K, 
which can be chosen arbitrarily or else completely scaled 
out of the problem. Indeed, if"/^ has units of length, 
time and mass, and has units of density in the geo- 
metrical units. Thus, in the following, we only show non- 
dimensional quantities, which are rescaled by K. In other 
words, we adopt the units with c = G = K = 1. 

Following Butterworth & Ipser (1976), the line element 
is written as 

= _e2-df2 ^ 52g-2.^2 gj^2 ^^^^ _ ^^^)2 

+e^^-^-{dr^+r^de^), (17) 

where u, B, to, and C are field functions. The first three 
obey elliptic-type equations in axial symmetry, and the 
last one an ordinary differential equation. These equa- 
tions are solved using the same method as that described 
by Shibata & Sasaki (1998). 

The total baryon rest-mass M*, Komar mass (gravita- 
tional mass) M, proper mass Mp, Komar angular momen- 
tum J, rotational kinetic energy T, and gravitational po- 
tential energy W are defined from the energy momentum 
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tensor or matter variables as 

M*=27r j pu*Be2^-2Vdrd(cos6l), (18) 

M = 27r J (-2T/ + T/)Be2^-2^r2rfrrf(cos6i), (19) 

Mp = 2n J pu\l + e)Be^^-^''r^drd{cos 0), (20) 

J = 2tt j phu*u^Be^'^-'^''r^drd{cose), (21) 

T=^jn, (22) 
W = Mp - M + T, (23) 

where W > 0. From these quantities, the well-known 
nondimensional parameters are defined as /3 = T/W and 
q = J/AP. The ADM mass Madm is defined from the 
asymptotic behavior of as 

Madm = - lim vr. (24) 

r— >oo 

For the stationary spacetime, Madm is equal to M (Beig 
1978). Thus, accuracy of the numerical solutions can be 
measured checking the relation M = Madm- 

In addition to these quantities, we often refer to the 
central density, pc, which is used to specify a rotating star 
for a given set of fi and M, and to the equatorial circum- 
ferential radius, R, by which a compactness parameter is 
defined as M/R. 

3.2. Analysis for secular st,ability 

The secular stability for rigidly rotating stars against 
quasi-radial oscillations can be determined by a turning- 
point method as established by Friedman et al. (1988) and 
subsequently used by Cook et al. (1992, 1994). Accord- 
ing to the turning-point theorem, a change of the sign of 
dM/dpc along a curve of a constant value of J indicates 
the change of a secular stability. Thus, in the numerical 
computation, curves of constant values of J are computed 
and plotted in the plane composed of M and pc for deter- 
mining the turning points. 

To determine the region for the stable stars, in the 
present case, one should compute dM j dpc along curves 
of a constant value of J. For {dM/dpc)j > (< 0), the 
stars are stable (unstable) against gravitational collapse. 
From this fact, one can distinguish the stable region from 
the unstable one. 

3.3. Numerical results 

Numerical computation was carried out using the spher- 
ical polar coordinates {r,6,(fi). For a solution of elliptic- 
type equations in axial symmetry, a uniform grid for r 
and cos6' is adopted with the typical grid size {Nr, Ng) = 
(1000, 160). To investigate convergence for the location of 
the turning points, the computations were also performed 
with Nr = 500, 750, and 1500 and Ne = 100. We found 
that the location depends very weakly on Ng as long as 
it is larger than 100. On the other hand, the convergence 
of the location is slow with increasing Nj. (cf. Figure 1). 
However, with Nr ^ 1000, the convergence is good, except 
for the region near mass-shedding limits. 



In the typic;al computations, the equatorial radius of a 
star is covered by 0.6A^,- grid points. In our method, the 
outer boundaries are located at a finite radius, and thus, 
the outer boundary conditions are approximate. However, 
the stellar equatorial radius is larger than several hundred 
M for stable rotating stars with n > 3 (cf. Table 2), imply- 
ing that the radius of the outer boundaries is larger than 
~ lOOOM. In this setting, even the approximate condition 
is very accurate. To check the magnitude of the numerical 
error associated with this approximation, we performed 
several computations changing the location of the outer 
boundaries with a fixed grid spacing, and indeed found 
that the errors in mass, density, and angular momentum 
are much smaller than 1%. 

Figure 1 is a summary of the numerical results for 
r = 1.328, 1.329, 1.330, 1.332, 1.333, and 4/3. In each 
panel, the dashed curve denotes the sequence of rotating 
stars at mass-shedding limits. Namely, there is no rigidly 
rotating star in the right-hand side of these curves. The 
solid curve denotes the sequences of constant values of J. 
The thick solid curve, the crosses, the solid squares, and 
the open circles with the dotted curve are sequences of 
the marginally stable stars with Nr = 1000, 500, 750, and 
1500, respectively. It is found that there is a critical value 
of J = Jmim where Jmin denotes a constant, below which 
all the stars are unstable. The values of Jmin are deter- 
mined at the intersection between the curves for the se- 
quence of the mass-shedding limit and for the sequence of 
the turning points (cf. Table 2). In the following, we refer 
to such a point as an intersection point. 

Here, we write the density of the marginally stable stars 
as pt{J) for J > Jmin- Prom the fact that all the nonro- 
tating stars are unstable for F < 4/3, we can determine 
that with p < pt for a given value of J > Jmin, the stars 
are unstable and otherwise, they are stable. Therefore, 
the rotating stars located between the sequences of the 
mass-shedding limits and of the turning points are stable 
against quasi-radial oscillations for F < 4/3. For F = 4/3, 
stars with p > pt for a given value of J are unstable. Thus, 
the stars located in the higher-mass side of the sequence 
of the turning points are stable. 

The long-dashed curves denote the relation of equation 
(8) with K = 1; the sequence of the marginally stable 
stars derived by the semi-analytic calculation in §2. It 
is found that the long-dashed curves agree approximately 
with the numerical sequences of the turning points for 
F > 1.329. The error in mass for a given value of the 
central density is typically 0.1%. This indicates that the 
semi-analytic calculation gives a good approximate solu- 
tion for the marginally stable stars. 

For F = 1.328, the long-dashed curve is located in the 
right-hand side (the higher-density side) of the dashed 
curve, where there are no rigidly rotating stars. Since in 
the semi-analytic calculation the configuration of the stars 
cannot be determined, unrealistic solutions for marginally 
stable stars may be derived. Thus, the long-dashed curve 
denotes the unrealistic sequence of the marginally sta- 
ble stars and indicates that there are no stable stars for 
F = 1.328. 

The intersection point denotes the most compact 
marginally stable star for a given equation of state. Unfor- 
tunately, it is not easy to identify the turning points accu- 
rately near the mass-shedding limits. This fact can be also 
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recognized from the fact that convergence of the numerical 
solutions for the marginally stable stars with increasing Nj. 
is achieved more slowly near the mass-shedding sequence 
than for other solutions. 

Figure 1 shows that the numerical sequences of the turn- 
ing points appear to bend sharply near the mass-shedding 
limits irrespective of the value of F. However, such ten- 
dency seems to be spurious because with increasing Nr, the 
point of the bending tends to approach the mass-shedding 
sequences and eventually disappear for Nr oo. This 
illustrates that the intersection point is not accurately 
determined from the raw data sets for Nr < 1500. In 
the present work, for the determination of the intersec- 
tion point, we pick up the data sets of the turning points 
slightly far away from the mass-shedding sequences, and 
perform an extrapolation of the turning points to the mass- 
shedding limits. The several quantities for the intersection 
point determined by this method are listed in Table 2. We 
estimate that the errors of the values for M, J, pc, and f3 
are - 0.1%, 1%, 10%, and 0.1%, respectively. 

Note that our results for n = 3 slightly disagree with 
those by Baumgarte & Shapiro (1999). According to our 
results, the values of pc and q at the intersection point are 
w 1.1 X 10"^ and 0.91, while their results are w 0.7 x 10"^ 
and 0.97, respectively. This disagreement comes partly 
from numerical errors contained in our results of magni- 
tude ~ 10% for Pc and ~ 1% for q. However, even if this 
error is taken into account, the results of the two groups 
disagree. Currently, the reason is not very clear. Here, 
we point out a fact: In the limit Pc ^ 0, the mass of 
the marginally stable stars from our results converges to 
« 4.555, which agrees with the exact value, the value of 
the Newtonian polytrope with n = 3. On the other hand, 
the asymptotic value of their results slightly disagrees with 
the exact value in this limit. Thus, as far as the results for 
Pc — > are concerned, our numerical results appear to be 
more accurate. 

Figure 1 shows that for the smaller value of F, the al- 
lowed region for the stable stars becomes narrower. We 
find that for F < 1.328, there is no region for the stable 
stars and that for F = 1.329, the region for the stable stars 
is very narrow. From these results, we conclude that all 
the rigidly rotating stars are unstable for F 1.329. 

According to a semi-analytic estimate presented in §2, 
marginally stable stars satisfy equation (12). In this 
relation, the terms 3r - 4, 2(3F - 5)/3, and 2(3F - 

h)k^k2^ M"^/^ py^ are of nearly identical order as ~ 10~^ 
(except for the cases F w 4/3 in which 3r— 4 w 0). To con- 
firm the validity of this relation, we calculated the relative 
error defined as 



Q = 



2(3r - 5)/? 



3r - 4 - 2(3r - 5) ( /? - ^M2/3py3 



(25) 

Figure 2 shows the values of Q as a function of J for 
F = 1.329, 1.330, 1.332, and 1.333. It is found that \Q\ 
is less than ~ 2% irrespective of F. This indicates fair 
validity of the relation (12). 

For F — 1.329-1.332, the value of Q appears to asymp- 
totically approach a constant with increasing J (with de- 
creasing Pc). This is reasonable, because even in the limit 
Pc ^ 0, a marginally stable star is rotating, and hence, the 
semi-analytic calculation in which the stars are assumed to 



be spherical contains certain systematic error. Therefore, 
the asymptotic value of Q may be regarded as a typical 
magnitude of the systematic error associated with neglect 
of the nonspherical deformation in the semi-analytic cal- 
culation. 

Equation (9) implies that if qmin is larger than unity, the 
value of q for all the marginally stable stars is larger than 
unity. As shown above and in Table 2, q^i^ is larger than 
unity for 1.329 ^ F ^ 1.332, implying that the value of q 
is always so for the rigidly rotating and marginally stable 
stars. 

4. PREDICTING THE FINAL OUTCOME 

The marginally stable stars determined in §3 may be 
plausible approximate initial conditions for a stellar core 
collapse or a pair-unstable collapse. Here, we predict the 
outcome of the collapse paying particular attention to the 
black hole formation case under the assumptions that (i) 
the collapse proceeds in an axisymmetric manner, (ii) the 
viscous angular momentum transport during the collapse 
is negligible, and (iii) the pressure and heating effect never 
halt the collapse. Actually, in a stellar core collapse with 
mass larger than ^ 4OM0 (Fryer 1999) and in a pair- 
unstable collapse with mass larger than ^ 26OAf0 (Fryer 
et al. 2001), the pressure support and the nuclear burning 
energy are so small that the final outcome is likely to be 
a black hole. The numerical analysis is carried out in the 
same manner as that of Shibata & Shapiro (2002). 

Since viscosity is assumed to be negligible during the 
collapse, the specific angular momentum j of each fluid 
element is conserved in an axisymmetric system. Here, j 
is defined as 

j = hu^. (26) 

Next, we define rest-mass distribution m^{j) as a func- 
tion of j, which is the integrated baryon rest-mass of all 
fluid elements with speciflc angular momentum less than 
a given value jo: 



.{jo) = 2n [ , 



pw*Be2^-2''r2rfrrf(cos6l). 



(27) 



Let us assume that a seed black hole is formed during the 
collapse and consider the innermost stable circular orbit 
(ISCO) around the growing black hole at the center. If j 
of a fluid element is smaller than the value at the ISCO 
(.?isco)j the element will fall into the black hole eventually. 
Now the possibility exists that some fluid can be captured 
even for j > jiscOj if it is in a noncircular orbit. Ignoring 
these trajectories yields the minimum amount of mass that 
will fall into the black hole at each moment. The value of 
JiSCO changes as the black hole grows. If j/isco increases, 
additional mass will fall into the black hole. However, if 
Jisco decreases, ambient fluid will no longer be captured. 
This expectation suggests that when jisco reaches a max- 
imum value, the dynamical growth of the black hole will 
have already terminated (i.e., before the maximum, the 
point for j = jisco will be reached). As Shibata & Shapiro 
(2002) demonstrated, prediction of the growth of the mass 
and spin by this method is in good agreement with a nu- 
merical result in the collapse with n = 3, suggesting that 
this is also a good method for n ^ 3. 

To analyze the growth of the black hole mass, we gener- 
ate Figures 3 and 4. Figure 3 shows TO*(j)/M* as a func- 
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tion of j/M, for T = 1.329, 1.330, 1.332, and 1.333. Here, 
we choose the marginally stable rotating stars at mass- 
shedding limits. That is, we choose the most compact and 
most rapidly rotating, marginally stable stars. 

Each panel (i) in Figures 4 (a)-(d) denotes = 
J(j)/m*(j)^ as a function of m*(j)/M*. Here, J{j) de- 
notes the total angular momentum with the specific an- 
gular momentum less than a given value jo and is defined 
according to 



J j<3o 



'r'^drd{cose). (28) 



Now, J(j)/'m,^(j)'^ and ra^:{j) may be approximately re- 
garded as the instantaneous spin parameter and mass of 
a black hole, since the baryon rest-mass is nearly equal 
to the gravitational mass for the soft equations of state. 
Therefore, the solid curve in each panel (i) of Figure 4 
may be interpreted as an approximate evolutionary track 
for the angular momentum parameter of the growing black 
hole. It indicates that with increasing the black hole mass, 
the spin parameter also increases. 

If we assume that m^{j) and are the instantaneous 
mass and spin parameter of the growing black hole and 
that the spacetime can be approximated instantaneously 
by a Kerr metric, we can compute jisco (Bardeen et al. 
1972; chapter 12 of Shapiro & Teukolsky, 1983). In each 
panel (ii) of Figure 4, we show jisco['™*(.?)i 9*0)] as a 
function of m*(j). The maximum of jisco is reached at 
= 'Ti*crit) where m*crit/Af* ~ 0.57, 0.77, 0.89, and 
0.92 for r = 1.329, 1.330, 1.332, and 1.333, respectively 
(Fig. 4, circles). The value of jisco is slightly smaller 
than j at m^{j) = m^ciit- The point for j = jisco (Fig- 
4, triangles) is reached at a value of m^{j), which is only 
slightly smaller than m*crit- This suggests that the mass 
of the black hole will increase to ^ m^crit- However, af- 
ter the maximum of jisco is reached, iisco/'™*(i) steeply 
decreases above this mass fraction. Thus, once the black 
hole reaches the point of m*(j) = m^critj it will stop grow- 
ing dynamically. Panels (ii) of Figure 4 shows that at this 
stage, 9, « 0.94, 0.87, 0.77, and 0.73 for T = 1.329, 1.330, 
1.332, and 1.333, respectively. 

It is interesting to note that for F = 1.329-1.332, the 
total value of q for the system is larger than unity. How- 
ever, as far as the inner region of the collapsing stars is 
concerned, it is much smaller than unity and, hence, a 
seed black hole is likely to be formed. The large value 
of q > 1 will be reflected in the large mass fraction of 
formed disks surrounding the central black hole. For the 
cases with F = 1.329, 1.330, and 1.332, the fraction of 
the disk mass will be ^ 40%, 20%, and 10%, respectively. 
The present numerical analysis indicates that the global 
value of q is not always a good indicator for predicting 
the final outcome of stellar collapse with soft equations of 
state. However, the conclusion drawn here is based on sev- 
eral assumptions. To confirm it more strictly, fully general 
relativistic simulations are necessary. 

After the dynamical collapse, a system of a black hole 
and surrounding disks in a nearly quasicquilibriiim state 
will be formed. Subsequent evolution of the system is de- 
termined by the viscous time scale of the accretion disks. 
Since a large fraction of the total mass forms disks that 
eventually fall into the central black hole, the black hole 



will spin up to q 1 for F <^ 1.332. 

Next, we pay attention to the fate of the collapse for 
general marginally-stable stars. For the analysis, we only 
need to repeat the same procedure described above. Fortu- 
nately, without detailed numerical computations, it is pos- 
sible to approximately predict the outcome using a semi- 
analytic calculation as follows. 

In the semi-analytic method adopted in §2, the mass and 
angular momentum contained inside a cylinder of radius 
zu are written as 



M{w) = MM{y), 
J{ru) = JJ{y), 



(29) 
(30) 



where y = wS^i / Rs and denotes the Lane-Emden radial 
coordinate of the stellar surface (e.g., Shapiro & Teukol- 
sky 1983). M. and J are computed from the Lane-Emden 
functions as 

M{y) = ^ r 9^e{l-cost)d^, (31) 
^1 Jo 



where t = sin ^{mm[l,w^i/{RsS,)]), O'l denotes \d6/d£,\ at 
^ = ^1, and 

2 

6i"^^rfe (33) 



a^^e'i Jo 

Now we define the nondimensional angular momentum 
parameter estimated for the fluid elements inside the cylin- 
der of radius w as 



Jjw) J 
9M = T77=T2 =9^- 



M{tu) 



(34) 



Thus, a curve of qizu) is determined for a given value of q 
in this model. 

To verify that this semi-analytic estimate for q{w) is a 
good approximation, we plot (j(n7) as a function of A4 in 
each panel (i) of Figure 4 (dotted curves). To plot these 
curves, we choose the same value of q for each correspond- 
ing numerical model presented in Figure 4. Note that the 
specific angular momentum is written as ro^il in this semi- 
analytic model. Thus, M{w) may be regarded as the mass 
distribution as a function of the specific angular momen- 
tum. 

Comparison between the dotted and solid curves in Fig- 
ure 4 confirms that the numerical and semi-analytic re- 
sults are in good agreement. This indicates that the semi- 
analytic method adopted here is appropriate for an ap- 
proximately quantitative study. 

According to equation (34), the curves of the semi- 
analytic results shift upward with increasing q. For a suf- 
ficiently large value of q, the value of q{'cu = 0) exceeds 
unity and, therefore, q{vcj) > 1 for any value of ru. Such 
large values of q are achieved for marginally stable stars 
of small compactness. In this case, it is natural to expect 
that even a seed black hole will not be formed in the grav- 
itational collapse. Thus, the collapse of marginally stable 
stars with 1.329 < F < 4/3 and with a small compactness 
will not result directly in a black hole but in a disk or a 
torus. The critical value of q above which a black hole 
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will not be formed ((/bh) is ~ 2.5 irrespective of F (see the 
solid and dotted curves for q = 2.517 in Figure 4(a) for an 
example). 

Using equations (8), (9), and (13), the compactness of 
the marginally stable stars can be written as a function of 
Q- 

M ^ fc2(3 - n)a 

For q = 2.5, M/R^ = 2.61 x lO"-*, 2.02 x lO"", 8.18 x 10-^ 
and 2.06 x lO'^ for F = 1.329, 1.330, 1.332, and 1.333, re- 
spectively. In the iron core collapse of massive stars, the 
radius and mass before the collapse are ~ a few thousand 
kilometer and at most a few Mq, respectively (Umeda & 
Nomoto, unpublished). Thus, the compactness is ~ 10""^ 
and, hence, the value of q at the onset of the instability is 
likely to be smaller than 2.5 if the star is assumed to be 
rigidly rotating. On the other hand, in the pair-unstable 
collapse, the density and mass at the onset of the collapse 
are ~ 10^ g/cm^ and IOOMq (Bond et al. 1984). This 
implies that the compactness is ~ 2 x 10~^. At the onset 
of a pair-unstable collapse, the adiabatic index of the star 
decreases from ~ 4/3 to smaller value on a time scale much 
longer than the dynamical time scale (Bond et al. 1984). 
If the star is rotating sufficiently rapidly, the instability 
will set in when F decreases below ~ 1.329. For F = 1.329 
with M/Rs ^ 2 X lO"**, the value of q at the onset of the 
instability is larger than 2.5 if the star is rigidly rotating. 
This suggests that the collapse may not lead to a black 
hole directly in this case. 

A collapse with q > qbh will result in formation of a disk 
or a torus, which will be subsequently unstable against 
nonaxisymmetric deformation. After the nonaxisymmet- 
ric instabilities turn on, angular momentum will be trans- 
ported from the inner region to the outer region, decreas- 
ing the value of q around the inner region below unity and 
resulting in the formation of a seed black hole. 

To clarify whether the above scenario is correct, it is ob- 
viously necessary to perform a three-dimensional numer- 
ical simulation in general relativity (e.g., Shibata 1999; 
Shibata & Uryu 2000, 2002; Font et al. 2002; Shibata et 
al. 2003). During the collapse, the typical length scale 
may change by a factor of 10^ from Rg/M ~ 10* to 1. 
To follow the collapse by numerical simulation, very large 
computational resources will be necessary and, thus, the 
simulation for this phenomenon will be one of the compu- 
tational challenges in the field of numerical relativity. 

Finally, we note the following point. As shown above, 
q{'uj) near the rotational axis is smaller than the global 
value of by a factor of 2 in the rigidly rotating case. This 
is the reason that the relation g > 1 is not likely to be a 
good criterion for no black hole formation. However, in the 
differentially rotating case or for stiffer equations of state 
with F ^ 2, the ratio q/qivj = 0) becomes smaller than 
2. For q/q{w = 0) J; 1, the global value of q may be still 
an approximate indicator for predicting the outcome. Pre- 
vious axisymmetric simulations in general relativity have 
not been performed taking an equilibrium with F ^ 4/3 
in the rigidly rotating initial condition (Nakamura 1981; 
Stark & Piran 1985; Piran & Stark 1986; Nakamura et al. 
1987; Shibata 2000). This is the reason why the previous 
numerical works have suggested that the value of g is a 
good indicator for predicting the outcome. 



5. SUMMARY 

We have investigated the secular stability of rigidly ro- 
tating stars against a quasi-radial oscillation in general 
relativity. It is found that all the rigidly rotating stars 
with F < 1.328 are unstable against a quasi-radial oscil- 
lation. Stars with 1.329 < 4/3 that are unstable for 
the spherical case can be stabilized by the effect of rota- 
tion. Therefore, for rapidly rotating stars for which the 
adiabatic index gradually decreases from ~ 4/3 to smaller 
values, the instability against gravitational collapse will 
set in when the value of F decreases below « 1.329 (not 
4/3). 

The marginally stable stars for F = 1.329 4/3 have been 
determined numerically. It is found that for 1.329 < F < 
1.332, the nondimensional angular momentum parameter 
q for all rigidly rotating and marginally stable stars is 
larger than unity. The value of q for the marginally stable 
stars increases with decrease of the compactness. There- 
fore, for a sufficiently small compactness, the value of q 
exceeds unity even for 1.332 < F < 4/3. These results are 
in good agreement with those derived by a simple semi- 
analytic calculation presented in §2. 

The final outcomes of axisymmetric collapse of the 
rigidly rotating and marginally stable stars are predicted 
assuming that the pressure never halts the collapse and 
that the mass distribution as a function of the specific an- 
gular momentum is preserved. It is found that even for the 
case g > 1, a black hole will form as a result of the grav- 
itational collapse if the value of q is smaller than ~ 2.5. 
This is due to the fact that the angular momentum in 
the central region of the rigidly rotating stars is not large 
enough and that the effective value of the nondimensional 
angular momentum parameter q for the inner region is 
smaller than unity. The outcome of the collapse in such 
cases will be a rapidly rotating black hole surrounded by 
a massive disk. For q ^ 2.5, even the central region has 
angular momentum large enough to prevent collapsing to 
a black hole. In this case, the outcome will be a disk or a 
torus. It will be subsequently unstable against nonaxisym- 
metric deformation. Once the nonaxisymmetric structure 
is developed, the angular momentum will be transported 
outward by the nonaxisymmetric effects, decreasing the 
value of q in the inner region below unity and eventually 
forming a seed black hole. 

Although the conjecture mentioned above is reasonable, 
the final fate of the collapse of marginally stable stars with 
g ^ 1 can be determined only by fully general relativistic 
simulations. In particular for q ^ 2.5, three-dimensional 
simulations are necessary since the collapsing star is likely 
to become unstable against nonaxisymmetric deformation. 
Since the length scale will change by a factor of 10* during 
the collapse, huge computational resources will be neces- 
sary for the simulation. The simulation for this problem 
will be one of the computational challenges in the field of 
numerical relativity. 

Numerical computations were in part performed on the 
FACOM VPP5000 machine in the data processing cen- 
ter of National Astronomical Observatory of Japan. This 
work is in part supported by Japanese Monbu-Kagakusho 
Grants 14047207, 15037204, and 15740142. 
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Table 1 

Structure constants of spherical polytropes 



r 


ki 


k2 


fc4 


^5 


a 


1.328 


1.7877 


0.63546 


0.92247 


1.1837 


0.41330 


1.329 


1.7816 


0.63613 


0.92169 


1.1876 


0.41571 


1.330 


1.7756 


0.63680 


0.92091 


1.1914 


0.41810 


1.332 


1.7637 


0.63812 


0.91934 


1.1991 


0.42285 


1.333 


1.7578 


0.63878 


0.91856 


1.2028 


0.42522 


4/3 


1.7558 


0.63900 


0.91829 


1.2041 


0.42600 



Note. — Values of fcj and a are computed from the Lane-Emden function. 



Table 2 

Parameters of marginally stable stars at mass-shedding limits 



r 


M 


J 


q 


/3 


R/M 


Pc 


1.329 


5.435 


52.0 


1.76 


0.00847 


2.2e3 


1.3e-10 


1.330 


5.179 


37.0 


1.38 


0.00856 


1.4e3 


6.3e-10 


1.332 


4.779 


23.6 


1.03 


0.00873 


7.2c2 


4.5C-9 


1.333 


4.613 


20.0 


0.940 


0.00883 


6.0e2 


8.8e-9 


4/3 


4.560 


19.0 


0.914 


0.00886 


5.6e2 


l.le-8 



Note. — Adiabatic index, M, J, q = J/M'^, /S = T/\W\, R/M, and pc of marginally stable stars at mass-shedding limits are shown with the 
units c = G = K = 1. Note that M* is nearly equal to M in all the cases (the relative difference is smaller than 10~*). 
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Fig. 1. — Gravitational mass M as a function of pc for fixed values of J (thin solid curves), for sequences of rotating stars at mass-shedding 
limits (dashed curves), and for sequences of marginally stable stars (thicls solid curves) for F = 1.328, 1.329, 1.330, 1.332, 1.333, and 4/3 with 
Nr = 1000. Rigidly rotating stars are located on the left-hand side of the dashed curves, and thus, for F = 1.328, all the rigidly rotating stars 
are unstable. The long-dashed curves denote the sequences of the marginally stable stars derived by a semi-analytic calculation in §2. The 
crosses, filled squares, and the dotted curves with the open circles denote sequences of the marginally stable stars computed with Nj. = 500, 
750, and 1500, respectively. 
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Fig. 2. — Relative error Q for marginally stable stars defined in equation (25) as a function of J for F = 1.329, 1.330, 1.332, and 1.333. 
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Fig. 3.- Fraction of the rest mass m*(j)/Mt as a function of j/Mt for the marginally stable stars at mass-shedding limits with T = 1.329, 
1.330, 1.332, and 1.333. (J = 52.0, 37.0, 23.6, and 20.0, respectively; see Figure 1.) 
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Fig. 4. — 5* = J(j)/m, (jr)^ and j'lsco /Af, as a function of m, {j)/Mr. for rigidly rotating stars at inass-shcdding limits near the marginally 
stable criterion, (a) T = 1.329 and J = 52 (q = 1.76), (b) F = 1.330 and J = 37.0 (g = 1.38), (c) T = 1.332 and J = 23.6 (q = 1.03), and (d) 
r = 1.333 and J = 20.0 {q = 0.94). The crosses denote the point at which Jisco reaches maximum, the filled triangles the point at which 
j = jiscO) s-iid the filled circles the point at which = 1. The dotted curves denote the semi-analytic results for the same value of q as in 
the numerical results. For F = 1.329, the numerical and semi-analytic results for q = 2.517, for which g(ro) > 1 and hence a black hole will 
not be formed, are also plotted. 



